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I. INTRODUCTION 



There is some interest in computing the potential between two point-like sources wit 



i the 



qq colour quantum numbers in the octet (adjoint) representation o 



flB 



the colour group 



a. 



However, conflicting results have been reported in the literature [6|, |7|] and the question seems 
to still need some clarification. 

In this paper we reanalyze the problem discussing, in the pure Yang-Mills (YM) theory, 
the structure of energy eigenstates in the presence of colour sources and derive explicit 
formulae expressing the singlet and adjoint potential in terms of the Feynman propagation 
kernel (sometimes also called the Schrodinger functional) computed in the temporal (A = 0) 



gauge. The theoretica 



developed in refs. [9KL1]. 



framework we shall use is the formulation of the temporal gauge 



In sect. [TT] we review the general formalism. In sect. IIHI we illustrate the constraints 
imposed by global colour rotations on the structure of the functional integral. The solution 
to the problem of characterizing energy eigenstates belonging to different colour represen- 
tations is given in sect. IIVI in the formal continuum theory. In sect. |V] we give explicit 
formulae for extracting the singlet and adjoint gg-potential from the knowledge of the Feyn- 
man kernel. In sect. |Vl]we illustrate how to reformulate the previous analysis in the lattice 
language, suitable for numerical simulations. A number of technical issues are discussed 
in the Appendices. In Appendix |A] we recall some relevant group theoretical formulae. In 
Appendix [B] we construct the projectors necessary to single out states with definite global 
colour transformation properties. In Appendix Owe recall the formulae necessary to decom- 
pose the [N c ] (g) [N c ] tensor product into the sum of singlet and adjoint representation. For 
completeness we give in Appendix [D] the expression that the Feynman kernel in the presence 
of gg-sources takes in the Coulomb gauge. 



II. GENERAL FORMALISM 

The Feynman kernel in the presence of colour point-like sources^ be 
colour representations takes in the temporal gauge the expression 



ul 



onging to arbitrary 



K(A 2 , {u 2 }- A l5 { Ul }; T)= I V^h)TZ{U h ) {u2}{ui} K{A u 2 \ A l5 T) , (1) 



I Go 
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where T>fi(h) is the Haar invariant measure over the gauge group, Go, of the topologically 
trivial, time-independent gauge transformation that tend to the identity at spatial infinity. 
In eq. flTJ we have defined 

rA(x,T 2 )=A 2 (x) 
/A(x,Ti)=Ai(x) 

where Sy^A, A = 0) is the YM action taken at A = and 



__ r A(x,T 2 )=A 2 (x) 

K(A 2 ,A l ;T) = / VAexp[-S YM (A,A = 0)} (2) 



VA= [] dA(x,t). (3) 

x,T 1 <t<T 2 



In eq. (JT]) we have set T = T 2 — Ti and made use of the definitions [12J 

Aj^(x) = UtA k U h (x) + iUfaUhW , k = 1, 2, 3 , (4) 

[4(x) = exp[a a nx)]G£ , (5) 

where the matrices A a (a = 1,2, . . . ,N% — 1) are the hermitian SU(iV c ) generators in the 
fundamental representation [N c ], normalized to tr [A a A b ] = 5 ab /2 and we have indicated with 
A a = — (A a )* the generators of the conjugate fundamental representation [N c ]. Finally, we 
have introduced the compact notation 

L 

K(U h ) {m}{ui} = [| [expiAj^x,-)]^ (6) 



"2 U 1 



to represent a set of L sources in the colour representations [j] localized at the points x^. 
In eq. (jUJ) we denoted with [A^] u j u j the SU(iV c ) generators in the representation [j]. The 
indices u{ and u\ run over values appropriate for the representation [j]. 

Since K(A2, {^2}; Ai, {ui}; T) is the matrix element of the (Euclidean) time translation 
operator exp(— HT), we can write its spectral decomposition in the form 

K(A 2 , {u 2 }; A 1? { Ul }; T) = Y, e"^ T ^(A 2 , {u 2 }) ^(A l5 { Ul }) , (7) 

k 

where the functional ip k {A; {u}) is the eigenstate of the Hamiltonian, %, corresponding to 
the eigenvalue E k 

H^ k {A,{u}) = E k ^ k {A,{u}). (8) 

The state functional ^>&(A; {«}) transforms, under time independent gauge transformations 
(U w G Go), as 

M^ Uw , {«}) = E K\U w ) {u}{u , } il) k {A, {«'}) . (9) 

{«'} 



Eq. (J9j is equivalent to the Gauss' law in the Hilbert space sector of states with L external 
colour sources, as specified above. This can be checked by taking the functional derivative 
of both sides of the equation with respect to w a (x) at w a (x) = 0. 

It was shown in [9J that the kernel K(A 2 , {u^}'-, A 1; {^i}; T) and the eigenfunctionals 
^(A; {«}) are gauge independent in the sense that, had we chosen to quantize the theory 
in a different canonical gauge, say F(A) = (typically F(A) = V • A), the corresponding 
kernel and the related eigenfunctionals in the F-gauge would coincide with those in the 
temporal one, when evaluated on a gauge field A satisfying -F(A) = 0. 

Since the quantity J2{ u } ^*(A, {u})(j)(A, {u}) is invariant under gauge transformations 
(see eq. (jl])), in order to define a scalar product we must make use of the Faddev-Popov 
(FP) procedure. This is done by choosing any spatial gauge-fixing condition, say F(A) = 0, 



and accordingly defining the scalar product by means of the formula 13] 

(V,$= /2?MA)X>*(A,{uMA,{u}), (10) 

M 

where 

Vfi F (A) = l[A F (A)5[F(A)]dA^) (11) 

is the gauge field integration measure. In eq. (ITT]) A^(A) is the FP determinant [14J defined 
by 

A F (A) / Vn(h)6[F(A u *)] = l. (12) 

JQo 
We recall that the value of scalar product in eq. ( ITU]) is independent on the gauge function 

F(A) chosen to define it. 

Denoting with dk the degeneracy of the energy level Ek, we can compute the complete 

trace of the Feynman propagation kernel ((7j) obtaining 

/ V^ F {A) J2 K(A, {u}; A, {u}; T) = £ 4e~^ T . (13) 

W * 

The gauge invariance of the l.h.s. of eq. (fi"3"]) implies the gauge invariance of the quantities 
dk and E^. 

III. GLOBAL COLOUR TRANSFORMATIONS 

In this section we discuss how global colour transformations are implemented in the 
temporal gauge when external sources are present. 



Let us denote by V a global colour rotation, i.e. a constant SU(iV c ) transformation. The 
global colour invariance of SVm(A, A q = 0) implies 

K{A v 2 ,A v l] T) = K{A 2l A l] T). (14) 

From the definition ([1]) we then have 

K(AV,{u 2 };AY,{ui};T) = 

= [ V^h)n(U h ) {u2}{ui} K(A v 2 u \AX-,T) = 

J Go 

= J V^h)K(U h ) {u2}{ui} K(A v 2 UhV \A l -T) = 

J Go 

= f V f i(h)n(V^U h V) {u2}{ui} K(A^,A l ;T) = 
J Go 

= K(V% 2}W2} f P/i(/i)W(Z7 fc ) M}{tll} ^(Af\A 1 ;T)W(10 { «iHt.i} = 

J Go 

= K(V% 2}W2 yK(A 2 , K}; A 1; K}; T)K(V) Wl}{ui} . (15) 



Notice that the functional change of variables VUhV^ — > Uh in the third equality of eq. ( IT5l) 
is allowed by the fact that VUhV^ is a (topologically trivial) gauge transformation that tends 
to unit at spatial infinity and hence belongs to Qq. 
Eq. ( fl5|) can be rewritten in the form 

TZ{V) {U2}W2} K{A V 2 , K}; A 1; { Ul }; T) = 

= K{A 2 , {u 2 } ; Af, K}; T)K(V) {U[}{U1} . (16) 

Eq. (fl~6l) allows us to define the operator U(V), implementing global colour transformations, 
which acts on wave functionals as 

[W(V>](A, {u}) = K(V) {u}{u , } ^(A v , {«'}) . (17) 

U(V) provides a unitary representation of the SU(iV c ) group in the Hilbert space of state 
functionals and commutes with the kernel K. 
The unitarity of U(V) follows from the chain 



(W^V^M) = &M(y)4>) = J V^ F {A)r{A,{u'})K{V) {ul]{u} <P{A v ,{u}) = 
= J Vf, F (A)[^(V) {u}{u ^(A v \{u'})r<t>(^{^}) ■ (18) 

In fact, eq. ( IT8l) shows that 

[U\V)i>](A, {u}) = n\V) {u}W ^(A v \ {«'}) , (19) 



which, together with eq. (117)1 . gives U(V^) = WiV) entailing the unitarity of U{V). 

To prove that the Feynman kernel commutes with U we take an arbitrary wave functional, 
ip(A, {«}), in the appropriate source sector (see eq. ©), and consider the chain of equalities 

U{V) J Vfji F {M)K(A 2 , {u 2 }; At, { Ul }; T)^(A 1; {u,}) = 

= K{V) {U2}Wi} fvn F {X x )K{Al, {u' 2 }- A l5 { Ul }- T)i/>(Ai, {m}) = 

= JVMM)K(A 2 , {u 2 }; AY\ «}; T)^(V) {<}{ui} ^(A l7 { Ul }) = 

= J V l i P {X 1 )K{A 2 , {u 2 }; A l7 {u'J; T)n{V) {<}{ui} ^{AX , {%}) = 

= j Vfi F (A 1 )K(A 2l {u 2 }- A l5 {«!>; T) [W(y)^](A 1; { Wl }) , (20) 

where the third equality is a consequence of the invariance of the measure (fill) under global 
colour rotations. Owing to the arbitrariness of ip(A, {«}), eq. (|2"U|) implies that U{V) com- 
mutes with the Feynman propagation kernel. 

Setting Ai = A 2 = A in eq. (fT5l) and integrating over A with the measure (fill) , we get 
the key formula 

K {U]{U , } (T) =Jvn F (A)K{A, {u}; A, {u'}; T) = J Vfi F (A)K(A v , {u}; A v , {u'}; T) = 
= K{V% }{V} K {V}{V , } (T)K(V) W}{U , } . (21) 

Eq. (121]) tells us that K{T) commutes with TZ(V), so that by Schur's lemma it is a multiple of 
the identity matrix within any irreducible colour source representation. If we are in a sector 
in which more than one colour source is present, 7Z(V) can be decomposed into the direct 
sum of irreducible representations of SU(iV c ), and K(T) itself is a direct sum of multiples of 
the unit matrix, one for each irreducible component. We will make use of this result below. 
We end this section with some considerations on what happens if the gauge integration 
over Q is extended to the set Qq of transformations that also includes global colour rotations. 
In this case the colour averaged kernel 

K(A 2 , {u 2 }; Ai, {M; T) = [_ V^{h)K{U h ) {u2]{ui] K{A u 2 \ A i; T) (22) 

J Go 
satisfies the property 

K(A%, {u 2 }; Ai, {m}; T) = ft(0) WK} lT(A 2 , {<}; A 1; { Ul }; T) , (23) 

that implies for the averaged state functionals, ip(A, {u}), appearing in its spectral decom- 
position, the invariance property 

[W(V$](A, {u}) = K(V) {U}W} ^(A V , {«'}) = ^(A, {u}) . (24) 



Eq. (j2"3|) means that, for the averaged state functionals, a global colour rotation on the gauge 
field amounts to a rotation TZ(V) acting only on source indices, or in other words that any 
such state functional is a global colour singlet. 



IV. THE 



-SOURCE SYSTEM 



In this section we specialize the formulation of sects. [TT] and IIHI to the particular case of 
two sources with the q and q colour quantum numbers and show explicitly how to define 
and extract the gg-potential in the singlet and adjoint representation. 



The Feynman kernel in this particular source sector has the form [9M11] 



K(A 2 , s 2 , r 2 ; A x , si, r t ; T) = 
V/j(h)\exp[i\ a h a (x q )} 



Go 



Qu 



Vfx(h) exp [i\ a h a (x q )] 



S2S1 



S2S1 



exp [t\ a h a (^)]\ K{A 2 U \ A i; T) = 

J T2T\ 

exp HA /i (x ? )]l K{A 2 U \ Ai; T) . 



r\T2 



(25) 



In order to pick up the energies of the lowest lying states one has to study the large T 
behaviour of the expression ( |25l) while, at the same time, projecting out the desired colour 
structure. To this end one needs to classify the energy eigenfunctionals in terms of their 
global colour transformation properties. 



A. Classification of the energy eigenstates 

Since, according to eq. (120]) . the Feynman kernel commutes with U(V), the energy eigen- 
states are classified in terms of the irreducible representations of the colour group SU(iV c ). 

Every energy eigenstate in gg-sector is described by a wave functional which can be 
parametrized as 



V^(A; s, r) = [<P{A)I + A a a (A)] sr = 0(A)/ + A a a (A) = ^(A) 



(26) 



in terms of the 1 + (AT? — 1) functionals 0(A) and a (A). The parametrization given in 
eq. ( 126]) holds for any one of the components of an irreducible multiplet. The totality of the 
members of the multiplet can be reached through global colour transformations, according 
to the formula 



U(V)ip{A) = 4> V (A) = VtP(A v )V^ = <p(A v )I + VX a V^ a {A v ) . 



(27) 



As it can be seen from eq. ( 127)) . the global colour transformation of a state is composed of 
two different contributions: 

• the "colour-spin" contribution coming from the action of the matrix V on the source 
indices 

• and the "orbital" contribution coming from the transformation A — » A v . 

As we said, any energy eigenstate must belong to an irreducible representation of the 
colour group, then if) (A) must span a unique irreducible representation for any value of A, 
when transformed as in eq. (1271) . Notice, however, that for A = the orbital contribution 
is suppressed and eq. (127)) becomes 

^(0) = V^(0)V^ = 0(0)/ + VX a V^ a (0) . (28) 

Eq. (128)) implies that the two terms in the r.h.s. cannot be simultaneously different from 
zero, otherwise "0^(0) would belong to the reducible I® [N% — 1] representation. This means 
that we have the following three alternatives 

1) a (O) = and 0(0) ^ 

2) 0(0) = and a (O) ^ (for some a) 

3) 0(0) = a (O) = 

in correspondence to three different types of irreducible representations that we are now 
going to discuss. 

1) - 2) When 0(0) ^ (or a (O) ^ 0), i/)(A) must belong to the "colour-spin" singlet (or 
adjoint) representation, because the representation it belongs to cannot change discontinu- 
ously under a smooth variation of the "parameters" A. 

3) For states such that 0(0) = a (O) = we can have more general assignments of global 
colour quantum numbers. 

Recalling that the adjoint representation R[ A< ^(V) is defined by 

VX a V^ = R^ d \V)X c , (29) 

with R^ Ad ^(y) real, eq. (127)) can be rewritten in the form 

^(A) = <P(A v )I + Ri A a d \V)M^ V )^- (30) 



The two A-functionals in the r.h.s. of eq. ( 130]) must transform under the colour group and 
therefore will display "magnetic" quantum numbers. Thus we will denote them as <f) m (A), 
4> a k(A) respectively, with transformation properties under A — y A v given by 

<p m (^ v ) = R 1 ±>(V)MA) , (31) 

<t> ak (A v )=R [ ®(V) ( f )ak/ (A), (32) 

where R^ and Rffi are SU(iV c ) representations. Using eqs. ( 13T1) and ( 132]) . eq. ( 1301) becomes 

1^(A) = rt±,{V)ct> m ,{A)I + R [ * d \V)Rf kl {V)<p ak ,{A) A c , (33) 

showing explicitly that the </> m 's transform according to representation [a], while the 4> a kS 
belong to the direct product [f3] <g) [N% — 1]. 

Due to the trace-orthogonality of the identity and the A- matrices, the two terms in eq. (|33|) 
must separately belong to irreducible representations. 

We conjecture that the two terms in eq. (133]) correspond to two different types of energy 
eigenfunctionals, although we cannot exclude that, when the two irreducible representations 
are equivalent, both terms may simultaneously contribute to a given eigenfunctional. 

It follows from this analysis that the [a] representation must be irreducible and that, 
owing to the colour invariance of the dynamics, the "wave function" $ ak will be such to 
single out an irreducible representation from the tensor product \J3] <g) [N% — 1]. 

In summary, depending on the structure of the source indices and the colour represen- 
tation to which the gluon wave function belongs, we can classify the energy eigenstates in 
four classes. 



Colour-spin singlet, orbital singlet states 

ds] 

J [S] 



Wfj(A)=0(A)/ (34) 



<P(A V ) = 0(A) 
with 0(0) ^ 0. 
Colour-spin adjoint, orbital singlet states 



^(A) = A a o (A) (35) 

6JA V ) = 6 J A) , 



with 4> a (0) 7^ for some values of a. 



Colour-spin singlet, orbital [a] states with singlet source indices, and belonging to the 
irreducible colour representation [a] (see eq. (jHTl) ) 

^ 1 [ 5](A) = [al m (A)J (36) 

^ a] m (A v ) = R [ lAV)<P [a] m '(A), 

with M m (O) = 0. 

Colour-spin adjoint, orbital [/3] states giving rise to an irreducible colour representation 
[a] e IP] <g> [N 2 C - 1] (see eq. PJ) 

^Li(A) = A> afe (A) (37) 



J a k [ 



A v ) = Rtl,{V)<l> ak ,{A) 



with a fc(O) = 0. We repeat that in eq. (1371) a fc(A) is a tensor in the [/3] <g> [iV c 2 — 1] 
space living in the [a] -invariant subspace with values of a's and fc's constrained so as 
to yield the desired value of the index m. 

We remark that, since A v = A for V G Zjv c , where Zn c is the center of SU(iV c ), the repre- 
sentation [a] in eqs. fl36l) and fl371) is actually bound to be a representation of SU(iV c )/.Zjv c . 
The above classification is confirmed by perturbation theory [lOj [15| [16 1. States like 
those in eqs. (|3"I|) and (|35|) correspond, at zero order in the coupling constant g, to wave 
functionals which do not contain any gluon component, and are therefore non vanishing 
at A = 0. Excited states consist of states with gluons added to external sources. The 
presence of gluons allows the energy levels to reach any zero iV c -ality (triality for N c = 3) 
representation. States of this kind, however, contain some power of the gluon fields, so 
that their wave functional vanishes at A = 0. As a working hypothesis, to be checked by 
non-perturbative numerical simulations, we conjecture that similar structures also describe 
the non-perturbative dynamics of the external source sectors. 

B. The structure of the A-traced propagation kernel 

In this section we determine the way in which the four classes of colour multiplets, 
discussed in sect. IIVAI contribute to the Feynman propagation kernel traced over the field 

10 



boundary values 

K S2r2;siri (T)=jv^ F (A)K(A,s 2 ,r 2 ;A,s h r 1 ;T). (38) 

We remark that the gauge fixing in the A integration in eq. (1381) is needed in order to make 
all the matrix elements of K S2r2 - Siri (T) finite. In fact, although K(A, S2,r 2 ; A, si,ri; T) is 
invariant only under gauge transformations which are equal to the identity at the location 
of the colour sources, the colour-traced kernel J2 sr -^(A, s, r; A, s, r; T) is invariant under all 
the time independent gauge transformations belonging to Q , and therefore J2 sr K sr]Sr {T) 
would be infinite in the absence of gauge fixing. 

The contribution of the states given in eqs. (|34|) . (1351) . (1361) and (1371) to the partial 
trace (138]) . can be computed from the expressions of the partial traces of the corresponding 
projectors given in Appendix [B] We find 

• Colour-spin singlet, orbital singlet states of the form ([3] 



^ 2 ;^(T)|^ = ^^e- slS]r . (39) 



Colour-spin adjoint, orbital singlet states of the form ( )35l) 



K S2rr , Siri (T)Q d] = 2j2\ C s 2 r 2 K lSl e~ ElAd]T • (40) 



Colour-spin singlet, orbital [a] states of the form (1361) 



if.,» n m|^^w™«- !l * (4i) 



Colour-spin adjoint, orbital [(3} states of the form ( 13T|) with [a] G [f3] <S> [N% — 1] 

Ks 2 r 2 ;s 1 r 1 \T) = 2— ^ 2^ A S2 r 2 A nsiC . (42) 



C. The structure of the K(0, S2,V2',0, si,ri;T) kernel 

The previous considerations immediately imply that only the states in eqs. ( 1341) and ( )35l) 
contribute to i^(0, S2, r 2 ; 0, si, r\\ T), yielding terms with a tensor structure proportional to 
that of eqs. (1391) and ( HU1) . respectively. In formulae we get (see eqs. ( 1B2J) and ( 1B5I) ) 

i^(0, s 2 ,r 2 ;0, Si,n;T) = 

= l0(O)r%^e-^ + El0 a (O)| 2 EAU^ 1 e- £lAd]T + ... ■ (43) 

iV c a b 

11 



From the structure of states summarized in equations from (IMj) to (1371) and the results in 



perturbation theory [10|] 15( 16|, we conjecture that only the first two terms displayed in 
eq. f !43p are actually present. This statement can be tested in lattice numerical simulations. 
If we were only interested in the computation of the gg-potential in the vacuum, the 
relevant quantity to study would be K(0, S2,r 2 ; 0, S\,ri] T). Instead, in order to explore 
more complicated situations such as, for instance, the gg-potential at finite temperature, the 
relevant quantity to consider is K S2r2]Siri (T). 

V. EXTRACTING SINGLET AND ADJOINT POTENTIAL IN THE CONTIN- 
UUM 

To extract the lowest energy eigenvalues of the singlet and adjoint channel we have to 
project eq. (138]) (or eq. (J4"3"|) ) over the S r2S2 6 siri /N c and 2J2 a K 2 s 2 K iri tensor projector (see 
Appendix [U]), respectively. 

Referring to the partially traced kernel (138]) and using eqs. (139]) and (HT]) . we obtain at 
large times for the "colour-spin singlet" channels (see eqs. (13~4"]) and ( 131?)) ) 

E ^ r2S2 5 Siri K S2r2 . tSiri (T) T 4°° e- ElS]T + ... + D [a] e~ Ela]T + ..., (44) 

S2V2S1T1 -'"C 

where the dots stand for exponentially suppressed terms and, as discussed before, D[ a -\ is 
the dimension of the colour representation [a] to which the energy eigenstate belongs. 

Similarly, using eqs. ( 140]) and ( 142]) . for the "colour-spin adjoint" channels (see eqs. ( 135]) 
and (13TJ) ) we find 



E aEW^mnm^^-lle-^ + .-. + V^^... • (45) 



In perturbation theory the "singlet and adjoint gg-potentials" defined in eqs (144]) and (145]) 
contain divergent self-energy contributions, independent of the relative position of the 
sources. In perturbation theory these contributions can be renormalized away dividing the 
Feynman kernel (13"8]) by the product of the two traces (characters) of the kernels with the 
insertion of a single q or q source. This can be done by computing, at large T, the quantity 

T (T) = ^r 2 ;s 1 r 1 (T) 

12 



where 



KJT)= f Vn{h)tr[e iXaha ^] f Vfi F (A)K(A Uh , A; T) , (47) 

JQo J 

Kg(T) = f Vfi(h) tr [e^ 4 *'^] J V^ F (A)K(A U \ A; T) . (48) 



We stress again that, had we included in the gauge average also the integration over global 
colour rotations, the states in the spectral decomposition of the averaged Feynman kernel 
corresponding to colour non-singlet states would be missing (see eq. ( 124]) ). Consequently, in 
the large T limit, we might not be anymore in position to reach the desired lowest energy 
eigenvalues. Another unwanted consequence is that colour averaging makes the Feynman 
kernel contribution of a state vanish when the tensor product of source representations 
does not contain the singlet. This is precisely what would happen to the factors in the 
denominator of eq. (145]) . 

In the non-perturbative regime, where colour is supposed to be confined, the quantities 
defined in eqs. (147]) and (148"]) could be zero anyway. So the need and the way in which 
this particular renormalization step should be carried out must be the object of numerical 
investigations. 

VI. EXTRACTING THE ^-POTENTIAL FROM LATTICE SIMULATIONS 

In this section we discuss how to compute, in lattice simulations, the singlet and adjoint 
gg-potentials. Although it is not difficult to translate in lattice language the temporal gauge 
fixing procedure presented in the previous sections and provide a discretized version of 
eqs. (138]) or (143]) . it turns out that these formulae are not well suited to the structure of 
practical numerical simulations. 

The reason is that it is not easy to implement on the lattice the condition Uh £ Go (i- e - 
Uh(x) — > I) on the boundary gauge integration in eq. (Op). The limitation Uh £ Go is, 
however, absolutely crucial, as we have discussed, in order to avoid the cancellation from 
the colour averaged kernel of all the global colour non-singlet states (if not dynamically 
confined) . 

The way out to this difficulty is to isolate in the Feynman propagation kernel on the 
lattice the contribution of states belonging to given representations, averaging with the 

13 



SU(A r c ) characters. This idea is based on the formula 
£ / Vfx F (A)K(A v , s, r; A, s, r; T) = 

sr ** 

= £ / ZM A) £ # (A", S , r) (# (A, S , r))*e- E[Qfc,T = 

= E/^(A) E ^i(^)* 1 (A, S ,r)(^(A, S) r))*e-^^ = 

sr k,m,m' 

= E tr [i?^(y)]e- £lrafe]T E^ x W(y)]e- eWT , (49) 

which shows that performing a global colour rotation on one of the boundary gauge field 
and taking the full trace of the kernel amounts to replacing the multiplicity factor appear- 
ing in eq. f lT3|) with the character of the irreducible representation to which each energy 
eigenfunctional belongs. 

A. Temporal lattice gauge fixing 

We start discussing how the temporal gauge fixing procedure can be implemented on 
the lattice. Link variables are generated by some Monte Carlo algorithm with weight 
given by (minus) the exponent of the (gauge unfixed) plaquette action, assuming periodic 
boundary conditions in the four dimensions. Lattice points are denoted by {n Q , n}, with 
n = 0,1, . . . , N T , n k = 1, . . . , N L , k = 1, 2, 3. 

On each configuration the temporal gauge fixing is carried out subjecting the gauge links 
to the gauge transformation 

ft(n, t) = C/ t (n,t-l;n,t)...C/ t (n,l;n, 2)£/ t (n, 0; n, 1) , t = l,...,N T . (50) 

ft(n,0) = I, (51) 

which sets Uo(n, 0; n, t) — I at t — 1, . . . , Nt — 1. Under the gauge transformation (1501 the 
spatial links transform as 

U k (n, t;n + k,t) -> U' k (n, t;n + k,t)= fi f (n, t)U k (n, t; n + k, *)ft(n + k, t) , (52) 

t = l,...,N T . 

Notice that, owing to eq. (f5"Tj) . the spatial links at t = are untouched. 

The time ordered product of the temporal links (open Polyakov lines) that represent (in 
the unfixed theory, see also eq. fIDlj) of Appendix [D]) the external q and q colour sources at 
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n q and rig, become 



N T 

]JU (n q ,t- l;n g ,t) = 
t=i 

N T 

= Y[Q(n q ,t-l)U^n q ,t-l;n q ,t)n\n q ,t) = n\n q ,N T ) = n\n q ), (53) 

t=i 

[[[Uo(n q ,t-l;n q ,t)\ = 



t=i 



t 



= [ J] fi(n,-, t - l)C/o( n $, * - 1; %, ^ f K £)] = Q(n q , N T ) = Q{n q ) . (54) 

t=i 

Configuration averages performed at this stage will inevitably also include an average over 
global colour rotations. In the next subsection we describe a method based on character 
projections which exploits this colour average to select the contribution to the Feynman 
propagation kernel coming from states belonging to a given global colour representation. 

B. Character projection 

As discussed above, we can easily compute on the lattice the colour averaged kernel (see 
eq. f )25p ). which in continuum notation reads 

K S2> r 2 -s uri (T) = [t>Hf{A) I T>p(h) \e ixaha ^} [ e ^°M* K(A 2 U \ A l5 T) . (55) 

J Jg L J S2Si L -lr 2 ri 

Within the (^-integration we can factorize the integration over Qq times the integration over 
the SU(iV c ) group, getting 

K S2 , r2 . Sl , n (T)= [ VVV S2S3 V r * 2r fv^ F (A)K(A v ,s 3 ,r 3 ;A,s 1 ,r 1 ;T), (56) 

where K(A, s 3 , r 3 ; A, s\, r±; T) is the original Feynman kernel in the temporal gauge defined 
in eq. ( 1251) . 

We will now show that the information about gg-potentials can be extracted computing 
the character-weighted kernel 

K b] S2 ,, 2;si , ri (T) = / VV(x b] (V))*V S2S3 V; 2 r 3 fv^ F (A)K{A v ,s 3 ,r 3 ;A,s 1 ,r 1 ;T), (57) 

where x^(^0 is the character of the representation [7]. In the following we will also be 
interested in the total trace 

K h \T)^Y,K h] s , r] sAT), (58) 

rs 
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for which, in virtue of its Qq invariance, we have the further benefit that there is no need for 
a lattice gauge fixing while computing the boundary A-integration. 

The proof of these facts is based on the orthogonality property of characters 

LvV( X W(V)y X W(V) = 6 [am , (59) 

•J bU (iV c ) 

where x^(^) an d X^OO are ^ e characters of the irreducible representations [a] and [/3]. 
To proceed we need to study the contribution to K^ S2tr2;sun (T) and its total trace of 
the four types of energy eigenstates (1M1) - (!371) discussed in sect. HVl 

• From (f34l) the colour-spin singlet, orbital singlet state contribution to K S2 ,r2;si,n(T) 



is seen to be 



7-1-1 (rr \ I 5 ! _ f VAifi eWt / fin \ 

K S2,r 2 ;s 1 ,r 1 { 1 ) [s] ~ d h],[S] jj e . (60) 



■tM 



Thus for the fully traced kernel we get 

I7\i\ 



~ -nQ = 5 hUS] e- ElS]T . (61) 



From ( 1B6I) we find that the contribution of colour-spin adjoint, orbital singlet states ([35 
to the character-weighted kernel is 



^wirfl = 2 / W(^))* E(^ c Vt) S2r ,(A: iri )*e^ [Ad]r = 

[Ad] JSV(N C ) c 

= 2/ Pn^^jj^^-^^A^JA^Je-^. (62) 

JSU(Nc) ac 

The r.h.s. of eq. f )62|) is different from zero only if [7] = [N% — 1]. Tracing with 
<5 S2Sl 5 r2ri the resulting trace over the A matrices provides the factor 6 ac /2 that brings 
in the character of the adjoint representation Da-Rl^ c_1 '(^0 = X^ Nc ~ 00- From the 
orthogonality of inequivalent characters we thus get for the fully traced kernel 



i^(T)|!l = / ^^(x W (^))*X [ ^- 11 (^)e- S[AdlT = S^^e-* 1 ^ . (63) 
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Recalling the form of the projector (1B9I) . we see that the contribution of the colour-spin 
singlet, orbital [a] states fl36l) to the character- weighted kernel is 



K 



h] 



■S2,r2;si,»'il-' ) 



[S] 



S2r2 riSl 



SU(JV C ) 



-DV{^\V)TY, / 2?/*,(A)^ m (A^)0M m (A)V 



_g[a]^ 



>* -£i^T 



•/su(iv c ) mm , J 

W( X H(y))* X H(V) = ^^5 [7] , [a] e- sHr , (64) 



S2r2 riSl 



N r 



SV(N C ) 



N r 



where we have used the orthogonality of the wave functionals 0^ m (A) for different 
Cartan indices. For the fully traced kernel we get 



K b \T) 



[S] 



niUcy 



-E^T 



(65) 



Finally using eq. (1B13I) and the orthogonality of <ft am (A) wave functionals for differ- 
ent Cartan indices, we find that the contribution of the colour-spin adjoint, orbital 
\P] states (eq. (1371)). composing into the irreducible colour representation [a] to the 
character-weighted kernel is 



*[7] 



[Ad] 



•" S2,r2;si,ri \-L ) 

= f VV( X W(V))*j:RW ck ,(V) fv^ F {A)^{A)<p ak {AY\l 2r2 {\\ iri re~ E[a]T ---: 

JS ' [J ( N <=> akk'bc J 

= 2 f VV(x h \V)yY,RtUv)K 2r2 (X b Siri re- Ela]T (66) 

•J oU ( I\r. ) J-L _ 



and is different from zero only if [7] = [a]. For the fully traced kernel we get 



K b] (T) 



SU(iV c ) 



VV(x»(V))'x [ " ] (V)e- E "' T = 6 b]M e-* 



.eI°>}t 



)[«. 



since J2bk Rbk-bk(V) = X (Y) is the character of the representation [ 



o 



(67) 
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C. Extracting gg-potential energies from lattice data 

To extract the interesting gg-potential energies from lattice data we propose to use the 
lattice version of (the large T-limit of) eqs. (I6TJ) and (|63|) for the singlet and adjoint repre- 
sentation, respectively. These formulae simply read 

N con f 



AS] 



K x (T) 



[ffl 1 "con; 

= T7 E tr t^K)] tr [fi(n*)] 



M T^oo e -£[S] T 



(68) 



— Y [Ad] I [Ad] 

>K X (T)\ 

V 7 I [S] lat 

i ^conf 



N n 



E (x [ ^(fi(noc))* tr[fit(n 9 )]tr[Q(n ? )] 



(f) 



+ 






+ 



(69) 



»con/ ^ =1 

where N con f (^> 1) is the number of gauge configurations that have been generated, and 
the dots represent corrections due to the finiteness of N con f as well as terms exponentially 
suppressed in the large T limit. In eq. ( 16"9~|) n^ denotes the lattice point at infinity, far 
from the location of the colour sources, and ^(noo) identifies the global colour integration, 
unavoidable in lattice simulations, mentioned at the end of sect. HTH (see also the step from 
eq. fl55l) to (156]) where V is what we are calling here f^rioo)). 

Naturally the character of the [S] representation is unit, while for the adjoint character 
we have x [M (V) = 2£ a tr [X a VX a V% 

For the q and q self-energy we find 

N t 

-i ly conj 



K m (T) 



K m (T) 



lat N conf 



E(x bl (n(*oo))*tr fit (n 



M 



+ 



-¥ e + 



-i ™ conf 



lat N 



conf i =1 



E (x M (n(Hoo))*tr fl(n 



(«) 



T— >oo —ijMy 



+ ... -4 e 



+ 



(70) 
(71) 



where, we recall, x M (^) = tr [V] and x [?] (K> = tr[T^]. 

VII. CONCLUSIONS AND OUTLOOK 

In this paper we have derived explicit expressions for the singlet and the octet (adjoint) 
potential between two static, point-like sources with the colour quantum numbers of a qq- 
pair (see eqs. (I6"8l) and (1691) ). They have a particularly transparent form in the temporal 
gauge. For completeness in Appendix [D] we provide the formulae valid in the Coulomb 
gauge. 
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We have discussed in sect. |TTT]the importance of limiting the gauge integration in eq. ([I]) to 
the gauge transformations that tend to the identity at spatial infinity, showing that otherwise 
only global colour singlet eigenfunctionals would contribute to the colour averaged Feynman 
kernel. 

In sect. I VI I we have shown how the temporal gauge fixing can be implemented in practical 
lattice simulations and how it is possible to extract the singlet and adjoint gg-potential energy 
from lattice simulations, weighting the fully traced Feynman kernel with the character of 
the representation one is interested in filtering out. 

Explicit numerical lattice simulations are under way [l7| to check to validity of the analysis 
presented in this paper and the viability of the formulae we have derived for extracting singlet 
and octet gg-potential. 
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the Galileo Galilei Institute for Theoretical Physics for the hospitality during the workshop 
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Appendix A: Completeness and Group Integration 

When the degeneracy of an energy level is due to the existence of a symmetry group, it 
is possible to describe the completeness sum over the degeneracy subspace in terms of the 
invariant group integration. The basic idea is easily explained as follows. 

Starting from an irreducible representation R(g) of a group G and a basis \i) of the 
representation space, we have the definition 

R(gM = E m\R(9)M = J2RM\J) • (Ai) 

3 J 

The completeness relation can be expressed as 

I = J2 \i)(i\ = D [R] I VgR{g)\k)(k\R){g) , (A2) 

■ JG 

with \k) an arbitrary representation state. In eq. flA2j) D\f>] is the dimension of the space 
spanned by the representation [R\. In fact, from (see ref. |l8| ) 

/ VgR jk (g)R lm (gy = ^ , (A3) 

Jg D[ R ] 
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we get 

D [R] f VgR{g)\k)(k\R\g) = D [R] £ / VgR fk (g)\j}(l\R rk (gy = £ \l)(l\ . (A4) 

JG j,l JG I 

The interest of eq. (1A2I) is that it allows to write the completeness relation by only knowing 
a single (arbitrary) representation state, \k). 

Appendix B: Colour Projectors 

We apply the formalism developed in Appendix |X] to write the colour projectors necessary 
to single out the states listed in eqs. (1M]1 . (1351) . (136|) and (l3~Tj) . 

• We start with the normalization condition for the colour-spin singlet, orbital singlet 
states of eq. (I3"4"|) which reads 

1 = JviM F (A)ti [|^J(A)V[fJ(A)] = A c |l? MF (A)|0(A)| 2 . (Bl) 

The projector over a colour-spin singlet, orbital singlet state is given by 

VH] = 0(A 2 )0*(A x ) <Wn-i , (B2) 

so that it contributes to the partially traced kernel (138|) with a term 

/©MA)|0(A)| 2 W riSl e-^ = ^l£i e -^T ; (B3) 

as it follows from eq. (IBlj) . 

• The normalization condition for the colour-spin adjoint, orbital singlet states of eq. ( )35|) 

is 

l = /^ F (A)tr[^(A)VfJ d] (A)] = 

= tr[A a A fe ]|PMA)0 a (A)0*(A) = i£/^(A)|0 a (A)| 2 . (B4) 

The projector over the colour-spin adjoint, orbital singlet states can be computed with the 
help of the group integration formula (1A3|) obtaining (we recall that R ab (V) is a real matrix) 



?>£] = (Nl - 1) f VV (V^ d] (A\)V^) S2r2 (V^ d] (A\)V% 

= (iV 2 - 1) / vv(vx a v^ S2r2 (vx b v% ri UA 2 )MMy = 

J SU(iVc) 

= (n 2 c -i)l vv R{i d \v)x: 2r2 R [ f b d \v)(\i i j*M^)MM) 

JSV NA 



SU(iV c ) 

J2M^2)MMTT,K 2 r 2 (K 1 r i r ■ (B5) 
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Using the normalization condition (1B4j) . we conclude that a colour-spin adjoint, orbital 
singlet state contributes to the partially traced kernel of eq. (12T1) with a term 

£ [vM^lM^EK^iK^Te-^^ = 2j2K 2 r 2 K lSl e~ ElAd]T . (B6) 

a c c 

As a check of eqs. (1B3I) and (1B6I) we can trace them over colour source indices obtaining 1 
for the singlet and (iV| — 1) for the adjoint state, in agreement with the fact that the trace 
of a projector is the dimension of the space over which it projects. 

• The normalization condition for the colour-spin singlet, orbital [a] states of eq ( 13T)j) is 
(no sum over m) 

1 = /2>MA) tr [^\JA)%\JA)} = N c J V^ F {A)\^ m {A)\ 2 . (B7) 

Eq. (1A2j) gives for the projector over the ip\% m multiplet the expression 

P$ = D [a] f VV^\ m {A v 2 )^ m {AXY5 S2r2 5 riSl (B8) 

for any fixed value of m. From the transformation properties ( )36|) and using eq. (1A3|) . we 



can further elaborate eq. (1B8J) with the result 



P$=D [a] f Wi?i(^)0 H m '(A 2 )( J R£L"(^)0 H m »(Ai))*<5 S2r2( 5 riSl = 

L > JSV{Nc) 

= Y, [Q] m (A 2 )0 [a] m (Ai)*^ 2r2( 5 riSl . (B9) 



Starting with eq. flB8p . we now compute the contribution of this state multiplet to the A 



partial trace of the kernel, finding 

D [a] f VV /DM A) ^ m {A v )^ m {A v y5 S2r2 5 riS1 e- E[a]T = 

— ^Ma a -e^t mirV\ 

- -j-j-o S2r2 d risl e . (BlUj 

To derive eq. fIBlOl) we have exploited the colour invariance of the measure T>np(A) and the 
normalization of the group measure, /su(iv c ) ^^ = 1- As before, the total trace gives the 
dimension of the representation space, -Dui- 

• In order to apply the same procedure to the colour-spin adjoint, orbital [a] states of 
eq. (1371) . we start from the normalization condition (no sum over m or k) 

l = /P/i F (A)tr[^(A)V[^( A )]=^E/^^(A)|0a fc (A)| 2 , (Bll) 
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and we then construct the state projector (no sum over k) 

r\% = D [a] f vv<p ak {A v 2 w bk {AX)R^\v)R [ ^ \v)K 2r2 {\ d SlTl y . (B12) 

This expression can be simplified by recalling the transformation properties <p a k{A v ) = 
i?J^/(y)0 a fe/(A) (see eq. fl37|) ) and the fact that the wave-functions (j) a k must be such to 
project out the irreducible representation [a] from the tensor product [f3] <g> [N% — 1]. We get 
from eq. flBT2|) 



vi a j d] = D [a] f pv r 4ii,(y)0 afc ^A 2 )(41k^)^(A 1 ))* J Rif ] (^)4f (v)K 2 r^ d s iri r 

J oU (-/VcJ 

= D [a] f W^KA 2 )0: fe KA 1 )i2g afe „(y)(Hg 6fe ,(y))*A^ 2 (A d Siri )* = 



SU(iV c ) 

E^(A 2 )0I fe (A 1 )5:A^ 2 (A 6 sm )*, (B13) 



S2''2 
a,fc 6 



where in the second line the matrix R^. a k'(V) is the irreducible [a] component in the tensor 
product R l k %(V)RM(V). 

Going back to eq. (1B12J) . we see that the contribution of this state multiplet to the A 
partial trace of the kernel is 



D [a] f VV fvM^)4> ak (A v )M^ V rR^ d] (V)R^ d] (V)K 2r2 (\ d Siri ye- Ela]T = 

•J oU \JS c) J 

= ^fHK^r^e-^ . (B14) 



One can again check that tracing also over source indices yields precisely Dr a i. 



Appendix C: Decomposing the [N c ] ® [N c ] tensor product 

The [N c ] <g> [N c ] tensor product representation acts on the space of complex matrices w sr , 
s, r = 1, . . . , N c , endowed with the scalar product 

(W 2 , Wx) = J2 w 2*sr W Ur = tr [w|wi] . (CI) 

sr 

From the normalization tr [A a A 6 ] = ^5 a b, we have the identity 

w-^trH = 2^A a tr[A a w], (C2) 



which implies 



w = ^ tr [w] + 2 J2 ^ ^ [A a u>] , (C3) 
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where In c is the N c x N c unit matrix. In components eq. ( 1C3|) reads 

^ 2 r 2 = Yl W Siri 5 S2Sl 5 r2ri = J2 W sm (-TT^ara^in + 2 I] A s 2 r 2 ^r lSl ) ■ ( C4 ) 

Eq. (IC4|) is the algebraic identity 



1 * r n 



The two tensors in the r.h.s. of eq. ( 1C5|) are precisely the two projectors onto the singlet and 
the adjoint representation, as it follows from the identities |l| 

2 11 

2 A 2 — 1 

iv c a iv c a 

Notice that within each irreducible representation space the two tensors (1C6I) and flC7[) act 
as unit operators of appropriate dimension. 



Appendix D: The qq- kernel in the Coulomb gauge 

For completeness we give in this Appendix the explicit expression of the Feynman kernel 
in the presence of external sources in the Coulomb gauge. The formulae collected in the pre- 
vious sections can be straightforwardly rewritten in the Coulomb gauge because, as proved 
in (9J, not only the energies (that are gauge invariant quantities) but also the eigenstates of 
the Hamiltonian are exactly the same in the Coulomb and the temporal gauge. In fact, in 
going from the Coulomb to the temporal gauge, the Feynman kernel stays the same as the 
two gauges are related by a change of variables in the functional integral 9]. 

Explicitly eq. (125]) takes the form 

K(A 2 ,s 2 ,r 2 ;A 1 ,s 1 ,r 1 ;T)= f Zty c (A) JJ ^o(x,£) (Dl) 

•>A(x,Ti)=Ai(x) x,Ti<4<T 2 



e 



where 



-s ym (a,a ) Texp r / w r)dr)l [Texp f - i / A (x ? , T)dr 



rir 2 



P/xc(A)= I] A c (A)5[VA]dA(x,t) (D2) 

x,Ti<t<T 2 

is the Coulomb gauge fixed integration measure and Ac (A) is the corresponding FP deter- 
minant. We note that, unlike the integration over the spatial components of the gauge field, 
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the integration over the temporal one in eq. (1D1|) is extended to include also the values at 
the boundary times. This is related to the fact that the temporal component of the gauge 
field, A , plays the role of a Lagrange multiplier which enforces the Gauss' law at any instant 



including the initial and final time 19J. This is why the exponentials in eq. (ID1I) do not 
wind up in the time direction. 

Setting now as in (1381) . A x = A 2 = A and integrating over A, again with the Coulomb 
integration measure, we conclude that the resulting quantity coincides with the expression 
displayed in the r.h.s. of eq. (|38|) from which singlet and adjoint potential can be extracted. 



[1 
[2 
[3] 
[4] 

[5] 

[6] 

[7 

[8 

[9 
[10 

[11 

[12 

[13 
[14 
[15 
[16 
[17 
[18 



L. S. Brown and W. I. Weisberger, Phys. Rev. D 20 (1979) 3239. 

L. D. McLerran and B. Svetitsky, Phys. Rev. D 24 (1981) 450. 

S. Nadkarni, Phys. Rev. D 34 (1986) 3904. 

G. T. Bodwin, E. Braaten and G P. Lepage, Phys. Rev. D 51 (1995) 1125 [Erratum-ibid. D 

55 (1997) 5853]. 

O. Kaczmarek, F. Karsch, P. Petreczky and F. Zantow, Phys. Lett. B 543 (2002) 41. 

O. Philipsen, Phys. Lett. B 535 (2002) 138. 

O. Jahn and O. Philipsen, Phys. Rev. D 70 (2004) 074504. 

E. V. Shuryak and I. Zahed, Phys. Rev. D 70 (2004) 054507. 

G C. Rossi and M. Testa, Nucl. Phys. B 163 (1980) 109. 

G C. Rossi and M. Testa, Nucl. Phys. B 176 (1980) 477. 

G C. Rossi and M. Testa, Nucl. Phys. B 237 (1984) 442. 

In our notation the gauge coupling is absorbed in a rescaling of the definition of the gauge 

field. 

G C. Rossi and M. Testa, Phys. Rev. D 29 (1984) 2997. 

L. D. Faddeev, V. N. Popov, Phys. Lett. B 25 (1967) 29. 

J. P. Leroy, J. Micheli and G C. Rossi, Nucl. Phys. B 232 (1984) 511. 

J. P. Leroy, J. Micheli, G C. Rossi and K. Yoshida, Z. Phys. C 48 (1990) 653. 

A. Guerrieri, S. Petrarca, M. Testa, in preparation. 

G Y. Lyubarskii, "The Application of Group Theory in Physics", Pergamon Press (New York, 

1960). 

24 



[19] Actually one of the two boundary Ao-integrations can be dropped. 



25 



